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Abstract. Non-local equations of motion contain an infinite number of derivatives and commonly 
appear in a number of string theory models. We review how these equations can be rewritten in the 
form of a diffusion-like equation with non-linear boundary conditions. Moreover, we show that this 
equation can be solved as an initial value problem once a set of non-trivial initial conditions that 
satisfy the boundary conditions is found. We find these initial conditions by looking at the linear 
approximation to the boundary conditions. We then numerically solve the diffusion-like equation, 
and hence the non-local equations, as an initial value problem for the full non-linear potential and 
subsequently identify the cases when inflation is attained. 
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INTRODUCTION 



Several string theory models such as the p-adic string and the cubic (super-)string field 
theory (CSSFT) flf lead to non-local field equations of motion that involve an infinite 
number of derivatives. These models give rise to interesting dynamics and have been 
considered in the context of cosmology [2], particularly inflation as well as bouncing 
universes. In general, we can write their action as follows: 
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where the d'Alembertian operator is 
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where the latter equality only applies for a Friedmann-Robertson-Walker universe. V( y) 
is the scalar potential and F(D) = — (1 +4% 2 <xn)e~ aD . The field equations are therefore 
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As it stands, the computation of the evolution of yf(t) seems to be a hopeless task as 
it involves an infinite number of derivatives. Indeed, for a general differential equation 
of the form 
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with finite n, we must specify initial conditions ty, \j/^ n l \ We may naively think 

that as n — > °o any solution is allowed, as we would be able to write y/(t) as a Taylor 
expansion whose coefficients can arbitrarily be chosen. While this is not the case, as 
the field equation itself imposes an infinite number of constraints on the allowed initial 
conditions [3], formulating an initial value problem (IVP) for the field's evolution still 
seems hopeless. This leads us to question whether we can reinterpret the problem in a 
helpful manner. 

Defining *P(f,r) = e~ raU yf(t), where r is an auxiliary variable, and differentiating 
with respect to r we obtain a diffusion-like equation [|4|] 
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Given that 1) = e aU y{t) and ¥(f,0) = yr(t), the equation of motion © leads to 
the non-trivial boundary conditions 
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In principle, the pde © can now be solved as an IVP by specifying an initial profile 
r) , which satisfies the boundary conditions ©, and evolving it forwards in time. For 
the diffusion-like pde © this is, however, a numerically unstable procedure. In this work 
we will solve © as an IVP, but the instability requires us to employ the regularisation 
methods described in Ref. |5|]. 

Finding the initial profile of *P(*i,r) is non-trivial. As an example let us take <§ 2 = 
and V(yA) = — yr/A, which is the case of the p-adic string, and also assume that 
^(t^r) = 0. Applying the d'Alembertian to the boundary condition and using the dif- 
fusion equation we find an infinite set of conditions that the initial profile *P(/i, r) must 
satisfy. 
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etc. This is equivalent to having to specify all the constraints on the allowed initial 
conditions that we mentioned before. It seems that the diffusion-like equation failed to 
bring us any advantage. In the next section, however, we are going to see that this infinite 
set of conditions is trivially satisfied when the equations are linearised which gives us 
hope to be able to go beyond the linear regime using the diffusion-like equation. 



LINEARISED EQUATIONS 

By linearising i// around a value A such that \ff(t) =A + (j>(t) and expanding the action 
© to second order, we obtain the new field equations 
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where U (0) = Uq — c0 — jm 2 (j) 2 and the coefficients are related to the original scalar 
potential V(\j/): Uq = V(A) +A 2 /2, c = -V'(A) -A, m 2 = -V"(A). As we have done 
before, we can now define r) = e~ mD 0(Y) from which follows the new diffusion- 
like equation 



with the linear boundary condition 
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Going back to our p-adic string example with ^ 2 = 0, V(y) = — \jr /4, and using A = 1 , 
applying □ successively to the boundary condition and using the diffusion-like equation 
we have the infinite set of conditions on the initial profile: 
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etc, and clearly, 4>(? i? r) = exp(r ln3) satisfies all of these. 

For a general model and potential, employing separation of variables, the solution for 
<!>(*, r) is 
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where (t) solves a local equation of motion 

□0 = -ft) 2 

where b = ceo 2 / (m 2 — 1) and CO 2 satisfies the characteristic equation: 
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This equation admits in general more than one root which can be complex, depending 
on the value of m 2 and whether E, 2 vanishes or not. General conclusions can be drawn 
from the reality of ft) 2 alone. In particular, if a root ft) 2 is real, so must be and we can 

define <p = 50 and b = Bb, where B 2 = y^Xa [l -4<^ 2 (1 +4aco 2 )] e a( ° 2 such that <p 
is a canonically normalised field satisfying the local equation □<£> = —co 2 (p — b. On the 
other hand, if ft) 2 is complex, its complex conjugate is also a root and we must search 
for solutions of the kind = Em 0m sucn that is real even though the various ra are 
not. Redefining the fields m in terms of <p m as before, and further expressing these 
in terms of the canonically normalised real fields % m and o m , (p m = (Xm + i^m) I >/2, 
<Pm = (Xm — iGm) / V% , we obtain their respective equations of motion 
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(17) 
(18) 



This system of equations describes what is called a "quintom". A set of two real fields 
where one is a standard field and the other is a ghost. 
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FIGURE 1. The p-adic string effective potential with p = 3. The solid and dashed lines represent the 
full potential and the quadratically expanded potential around the hilltop, respectively. 

THE /7-ADIC STRING 

We will now look at two examples of the evolution of y^Y) . We will start with the p- 
adic string for which we have y 4 = m A s p 2 /g 2 (p — 1), A = 1, V(v) = —y p+l /(p + 1), 
<^ 2 = 0, and a = ln(p)/m 2 , hence, the characteristic equation is simply e a(0 = m 2 , 
with m 2 = pA p ~ l and provided the solutions are (XG) 2 = lnm 2 ± 2nni, i.e., it 

admits one real solution and an infinite number of quintoms which imply that an infinite 
number of initial conditions must be set, in agreement with Ref. @]. In this work we 
will assume than only some of the degrees of freedom are dynamical and consider 
separately the evolution of ^(t, r) for individual values of n. The effective potential 
Veff(v) = f(V(\j/) + \j/ 2 /2) is illustrated in Fig.dJ 

In Fig. [2] we show the evolution of i// with initial conditions very near the hilltop. 
For the real field (n = 0), *P(?i,r) = 1 - £e aco R r and = with £ = 1(T 3 , i// 

decays to the minimum of the effective potential and oscillates . We also show the 
evolution of the Hubble rate and compute the number of e-folds of inflation, which 
can be seen to be sufficiently large to solve the horizon problem. For the first quintom 
(n = 1), *P(fi,r) = 1 — £e aa R r cos(aa)fr), and v j'(? i ,r) =0, we see that the field performs 
oscillations around the hilltop that grow in time and are propagated to the evolution of 
the Hubble rate. Eventually the amplitude of the Hubble rate is large enough to cause 
this to vanish and a recollapse occurs. We also show the evolution of the field had we 
considered only the quadratic approximation of the potential at the hilltop. 

THE CSSFT STRING 

In this section we will look at the dynamics of the field in the cubic superstring field 
theory. Here we have y 4 = m^/g 2 , X = -1, V(y) = ^ 2 = 0.96, a = 1 /4m 2 with 
characteristic equation e aa) (1 — 4^ 2 aco 2 ) = m 2 , and m 2 = pA p ~ l . In this case, if A ^ 0, 
the roots of the characteristic equation are of the form aoo 2 = 1/4^ 2 + LambertW(n,x) 




FIGURE 2. The left panel shows the evolution of the field \y, the Hubble rate and the number of e-folds 
of inflation in the p-adic string with p = 3 and n = 0. The right panel shows the same quantities for the first 
quintom, n = 1 . Solid and dashed lines represent the full non-linear and linear evolutions, respectively. 
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FIGURE 3. The effective potential of the CSSFT string. The left panel shows the regions where the 
dynamics is described as single real field (solid line) or as quintom (dashed lines), when n = 0,-1. The 
right panel shows the same regions for the remaining values of n. 



where x= — m 2 e~ 1//4 ^ 2 /4<^ 2 . If A = 0, then aco 2 = l/4£ 2 . InFig.[3]we show the effective 
potential of this model and indicate the values of ty/ where the roots are real and where 
they are complex (and therefore the dynamics is described by quintoms). It is clear that 
in this case the minimum of the potential only admits quintoms. 

The evolution of the field with n = is shown in Fig.|4]for initial conditions *P(/i, r) = 

£e aco R r and ^(ti,r) = 0. This corresponds to a real field evolving initially. We see that 
at late times, unlike in the p-adic string, the Hubble rate does go to zero indicating that a 
quintom has appeared. This is expected as only quintoms are expected at the minimum 
of the effective potential. For the n = 1 case, the evolution is also distinct from what 
happens in the p-adic string. Here the quintom decays towards the maximum of the 
potential and inflation lasts forever. 
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FIGURE 4. The left panel shows the evolution of the field \y, the Hubble rate and the number of e-folds 
of inflation in the CSSFT string with n = 0. The right panel shows the same quantities for n = 1. Solid 
and dashed lines represent the full non-linear and linear evolutions, respectively. 
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